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ABSTRACT 

The Hubble Tarantula Treasury Project (HTTP) is an ongoing panchromatic imaging survey of stellar pop- 
ulations in the Tarantula Nebula in the Large Magellanic Cloud that reaches into the sub-solar mass regime 
(< 0.5 Mq). HTTP utilizes the capability of HST to operate the Advanced Camera for Surveys (ACS) and the 
Wide Field Camera 3 (WFC3) in parallel to study this remarkable region in the near-ultraviolet, optical, and 
near-infrared spectral regions, including narrow band Ha images. The combination of all these bands provides 
a unique multi-band view. The resulting maps of the stellar content of the Tarantula Nebula within its main 
body provide the basis for investigations of star formation in an environment resembling the extreme conditions 
found in starburst galaxies and in the early Universe. Access to detailed properties of individual stars allows 
us to begin to reconstruct the evolution of the stellar skeleton of the Tarantula Nebula over space and time with 
parcsec-scale resolution. In this first paper we describe the observing strategy, the photometric techniques, and 
the upcoming data products from this survey and present preliminary results obtained from the analysis of the 
initial set of near-infrared observations. 

Subject headings: Magellanic Clouds — stars: imaging — galaxies: star clusters: individual — galaxies: star 
formation 



1. INTRODUCTION 

The Tarantula Nebula (a.k.a. 30 Doradus, hereafter 
'30 Dor") in the Large Magellanic Cloud (LMC) is one of 
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the most famous objects in astronomy. The first astronom- 
ical refere nces to the Tarantula Nebula are more than 150 
years old (Ide la Cailldl 17551: lHerschel|[T847l) . For more than 
a decade the Tarantula's ionizing cluster, R136, was thought 
to be the most massi ve known single star (250-1000 M0, 
iFeitzinger et al.|ll980l). Apart from a few d issonant voices 
(e.g. lFeastll9'5llWalbornll973l:lMoffatll982h . the true nature 
of R136 continued to elude the majority of the astronomical 
community until Weigelt & Baier ( 1985) resolved it into eight 
components by using holographic speckle interferometry. 

The Tarantula Nebula and R136 are now considered to be 
the closest (and only known) starburst in the Local Group. 
The size of the Tarantula Nebula {^ 200 pc in diameter) and 
"local" density of OB stars parallel those observed in systems 
characterized by very intense star formation, such as the star- 
burst knots observed in interacting galaxies in t he Local Uni- 
verse and young galaxies a t high redshift (z > 5. | Meurer et al.l 
[19971 [ShaplevetalJl2003l:lHeckmanetal.ll2004l) . Thus dis- 
secting the stellar populations and inferring an accurate de- 
scription of the anatomy of the Tarantula Nebula allows us to 
reconstruct for the first time the temporal and spatial evolution 
of a prototypical starburst on a sub-cluster scale. 

In this paper we present an introductory overview of the 
"Hubble Tarantula Treasury Project" (HTTP PI: Sabbi, GO- 
12939), an ongoing treasury program designed to observe the 
entire region with the Hubble Space Telescope (HST) in the 
near-UV (F275W and F336W), optical (F555W, F658N) and 
near-IR (Fl lOW and F160W). The program is built on an ex- 
isting HST monochromatic survey in the F775W filter (GO- 
12499, PI: Lennon), designed to measure proper motions of 
runaway candidates. 

Bias-subtracted and flat-fielded images, as well as all the 
raw data, are immediately available to the public and can be 
downloaded from the Mikulski Archive for Space Telescopes 
(MAST). Photometric catalogs and mosaicked images for all 
the filters will be distributed by our team to the community 
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Figure 1. The Tarantula Nebula as seen with HST (Credit NASA, ESA, Lennon; GO-12499). The image is the result of 15 HST pointings with the F775W filter, 
using WFC3 and ACS in parallel. The projected area in the sky is ~ 14' X 12'. The dynamic range of the image has been compressed to show both bright stars 
and faint filaments of gas. Colors are from ESO WEI exposures in the Ha (in red) and SlI (in blue) filters. The ionizing cluster R136 is located in the middle 
of the upper half of the image, in the brightest (white) part of the nebulosity. The small bright cluster to the right of Rl 36 is Hodge 301, and the loose stellar 
association in the middle of the lower half of the image is NGC 2060. Once observations are completed, HTTP will cover approximately the same region in six 
additional filters, namely F275W, F336W, F555W, E658N, El lOW and F160W. 
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in a timely fashion. The survey progress can be followed 



through the website http://30dor.stsci.edu 

This paper is organized as follows: we present the motiva- 
tions for the survey in Section|2] and the observing strategy in 
Section |3] As of the writing of this paper, we have collected 
half of the data in the near-IR (NIR). The data reduction and 
the analysis of this portion of the dataset are discussed in Sec- 
tion in while the color-magnitude and color-color diagrams 
(CMDs and CCDs respectively) are presented in Section |5] 
A discussion of the results and conclusions can be found in 
Section|6] 

2. TRACING THE EVOLUTION OF A STARBURST 
2.1. The Cosmological Context 

Starbursts are short-l ived p eriods of intense, massive star 
formation (Searle et aI1ll973l) that occur in compact regions 
(10 to 1000 pc) and dominate th e overall luminosity of 
their host galaxy jHeckmanl l2005h . In their Annual Re- 
view iKennicutt & Evans! (120121) reminds us that the star for- 
mation rate per unit area in a starburst is much larger than 
in disk-averaged star formation rate surface densities (^ 
0.1 M0yr"^kpc~^). Evidence for recent starburst activity is 
found in 6% of the galaxies in the Local Universe (z -C 1; 
iLeeet all [2007. 20o |) and in 15% of those at h igh redshift 
(z > 1: lo'ConnelllT2005l: iDouglas et al.1 120101) . indicating 
that starbursts are a relatively common phenomenon. The 
impact of a major starburst on the evolution of a galaxy is 
dramatic, because it will shape the galaxy dynamics, stel- 
lar evolution and chemical compositions in ways that are 
largely dictate d by the intensity and the duration of the burst- 
ing event (i.e.lDekel & S i lkl [19861: iMa c Low & Ferr£jal lT999t 
IKennicutt & Evansll2012ll 

Starbursts are studied on two complementary fronts. High- 
redshift surveys collect information over a range of long look- 
back times for a large sample of objects, but suffer from lim- 
ited or no spatial information and crude age diagnostics. Stud- 
ies in the Local Universe, on the other hand, are affected by 
the small size of the sample, but can probe discrete structures 
within the star-forming region using more precise population 
diagnostics (Hunter & Elme green 2006; Tolstoy, Hill. & Toil 
I2009t iMcOuinn et alJl2012l) . In both approaches the charac- 
terization of starburst properties is still rudimentary and some- 
times contradictory. 

Studies in individual nearby dwarf galaxies and in surveys 
at redshift higher than ~ 1.0 indicate that a starburst may last 
for 5-10 Myr, suggesting that the violent feedback from the 
massive stars formed dur ing the burst suppresses f or a while 



any further star fomition (iTosi et al.lll989l iFerguso n & Babul] 
I998l:lAnnibaU et alJl2003llStinson et alj2007l:lAnnibaU et a!] 



20091) . Other studies in the Local Universe, however ad- 
vocate for much longer star-forming episodes (> 100 Myr; 
iKonstantopoulos et al. 2009; McQuinn et al. 2010), although 
it is not clear if their approach allows for sufficiently high 
temporal resolution. Since the spectral energy distributions 
(SEDs) of galaxies at high redshift can be interpreted only 
by assuming a duration for the burst, understanding how and 
for how long star formation can be sustained during a star- 
burst event has profound implications on our picture for un- 
derstanding galaxies. 

2.2. The Tarantula Nebula: A Unique Opportunity to 
Decipher Starbursts 

The Tarantula Nebula is by far the most luminous 
and massive star-forming region in the Local Group 



( IKennicutt & Hodg3 II986I) . Covering an area of ^ 
40, 000 pc^, the Tarantula Nebula is the closest extragalac- 
tic giant Hll region and is comparable in size to the unre- 
solved luminous Hll complexes observe d in distant galaxies 
(lOey et al.ll2003l:lHunt & Hiras hita'2009). 

Wit h more than 800 spectroscopically confirmed OB 
stars (lEvans et al.l 1201 ll) . some of which are a mong the 
most massive candidates dCrowtheret al. I | 20I0l). 30 Do r 
is often described as a mini-starburst (ILeithererl 11998b . 
The UV flux coming from 30 Dor is such that it has 
been proposed as a s mall scale local a nalog to the 
Lyman-Break Galaxies ([Meurer et aP 119971: IShapley et al.l 
120031: [Heckman et al.l 120041). Furt h ermore several authors 
(e.g. lO'Connell & Manganol [19781: lO'Connellet aH [1991 
iBrandl et al.l 120041) have suggested that the luminous knots 
observed in bursting galaxies such as M82 are made of (mul- 
tiple) 30 Dor-like systems. 

By virtue of its location in the LM C (~ 50 kpc; 
iPanagia et anil991l:|Pretrzvnski et al.ll2013l) . with Hubble we 
can resolve the Tarant ula Nebula into sing le stars down to the 
sub-s olar mass regime (I Andersen et al.l20 09: De Marc hi et al. 
2011^. The low inclination angle (^ 30°; .Nikolaev et al. 



2004) limits the line of sight confusion, and the foreground 
reddening is low because of the high Galactic latitude. 

The Tarantula Nebula is a very d ynamic region bot h in 
terms of run away-star candid ates (e.g. lEvans et al.|[201 ll) and 
gas motion |Ch u et a nil994'). Since it is a multi-stage star- 
forming region (Walbom & Blades 1997), where loose asso- 
ciations and very dense star clusters of different ages coex- 
ist, the Tarantula Nebula is a noteworthy window in which to 
study the different modes of star formation and their recip- 
rocal interaction. In particular, observational evidence shows 
that over the last > 25 Myr the star-formation history of the 
region (Figure [T) has been complex : 



' IWalborn & BladesI (119971) identified multiple genera- 
tions of stars ranging from < 0.5 Myr to ^ 25 Myr 
in 1' - 2' from the core of 30 Dor, R136. 



• iGrebel & Chul ([2000l) derived an age of - 20 - 25 Myr 
for Hodge 301, a cluster ^ 3' to the northwest of Rl 36. 

• NICMOS observations at the "frontier" between 
Hodge 301 and R136 showed that this region is strewn 
with several massive O stars embedded in dense knots 
of dust, suggesting that stellar feedback from the two 
clust ers may be triggering new episodes of star forma- 



tion (iBrandner 



; tnggermg 
et al.ll2001l) 



• ~ 7' to the west of R 136 the super-bubble created by 
the -- 10 Myr old OB association LH99 is filled by the 
expanding supernova remnant N 157B, which cont ains 
the most energetic pulsar known (IChen et al.ll2006l) . 

• Deep HST optical and NIR observations suggest that 
the core of 30 Dor, R 136, is an interacting double clus- 
ter (iSabbi et al.ll2012l) . 

As these previous studies hint, the Tarantula Nebula region 
is complicated. While there have been attempts to study in 
detail some of its specific parts, much o f the focus with HST 
to date has bee n on i t s ce ntral cluste r (IWalborn et al.lll999l 
Brandner e tHl I200II: Imi born et all 120021: [Andersen et all 
2009i:iDe M'archi et al.l201 ItiSabbi et al.l2012l) . These studies 



have yielded significant results, nevertheless a comprehensive 
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Figure 2. Passband of the filters used in HTTP as a function of wavelength. 

understanding of the Tarantula Nebula in its entirety can only 
be reached by detailed studies that include other star clusters 
and surrounding field stellar populations. 

2.3. The HTTP Science Drivers 

The ongoing HST panchromatic survey of the Tarantula 
Nebula HTTP (PI Sabbi, GO-12939) is built on an exist- 
ing HST program (PI Lennon, GO- 12499) in filter F775W 
(- SDSS i-band). HTTP was awarded 60 orbits of HST time. 
Figure|2]shows the wavelength coverage of the survey, once it 
is combined with the GO- 12499 dataset. 

The decision to survey the region over such a wide wave- 
length baseline was driven by several factors. Adding an opti- 
cal filter to the F775W dataset allows us to reconstruct the star 
formation history of the entire region and infer the changes in 
star formation rate with time down to ^ 0.5 M©. We pre- 
ferred the ACS F555 W to the wider F606W filter, because the 
latter, being quite red, would have provided a limited color 
baseline relative to the F775W filter Furthermore the Ha line 
falls in the F606W bandpass and in Hll regions this can affect 
the accuracy of the photometry. 

Supplementing our dataset with two NIR filters allows us to 
extend our study to the dustiest regions, in order to better con- 
strain the properties of cool stars such as red supergiant, red 
giant (RGB) and asymptotic giant branch (AGB) stars. Even 
more important, only NIR observations can identify regions of 
ongoing star formation that are still embedded in their dusty 
cocoons (see Figure|3]l. 

Previous NICMOS observations have shown that star for- 
mation took place in the recent past at the gaseous int erface 
between R136 and Hodge 301 jBrandneret al.l 120011) . By 
extending this study to the entire Tarantula Nebula we can 
identify pockets of ongoing star formation far from the main 
clusters and test the constructive properties of s tellar feed- 
back predicted by the mod els (e.g. lElmegreen &~L ada 1977; 
iVanhala & Cameronlll998l) . We have selected the FllOW and 
F160W filters, since they provide the highest throughput with 
the IR channel. 

Observations in the near-UV (NUV) have been designed to 
infer the temperature of the hotter stars and, once combined 
with optical measurements, constrain the amount of extinc- 
tion (IZaritskv.l999;iRomaniello et al. 200 2). We have chosen 



the WFC3 filter F336W because of its high throughput and be- 
cause it is specifically designed to constrain the Balmer break. 
To better determine the temperature of the hotter stars and to 
provide the largest baseline possible to characterize the ex- 
tinction law we added the even bluer filter F275W. 

Finally we are observing the Tarantula Nebula with the 
narrow-band ACS filter F658N (corresponding to Ha). Young 
(< 5Myr) low-mass pre-main sequence (PMS) stars can be 
easily identified in the optical CMD as a population of faint 
red objects, and can be used to trace the extent of young 
star- forming regions (Sabbi et al. 2007; Cignoni et al. 201(3; 
iGouliermis et al.ir2012HSabbi et al.ll2012h . With the aging of 
a stellar population, however, the colors of the PMS stars be- 
come too blue to be distinguished from the low-mass MS stars 
of the field, an effect that becomes even more severe when ex- 
tinction is highly variable, as is the case for the Tarantula Neb- 
ula. In prev ious studies (IPanagi a et al. 2000; De March i et alJ 
I2OIOII2OIII) we have shown that a combination of broad- and 
narrow-band filters can be used to identify low-mass PMS 
stars with active mass accretion, and that this approach is very 
efficient in picking up older (> lOMyr) and more diffuse 
episodes of star formation. In addition, the same strategy will 
allow us to identify main-sequenc e and evolved stars with Ha 
excess, such as Be and B[e] stars (lGrebeil997h . 

In summary the high sensitivity, spatial resolution, and 
broad wavelength coverage of HTTP will allow us to recon- 
struct the Tarantula's star-formation history in space and time 
on a parsec scale. In particular we will be able to (1) charac- 
terize when and where star formation occurred; (2) establish 
the length and the strength of the star-forming episodes and 
their spatial scale; (3) depict the life cycle of star clusters; (4) 
portray the role of stellar feedback in shaping the evolution of 
a starburst; and (5) probe the universality of the stellar initial 
mass function. In addition, a Bayesian fitting of the stellar 
SEDs will allow us to accurately determine R\/, and therefore 
supply information about the size of the dust grains and their 
spatial distribution (K. Gordon et al. 2013, in preparation). 

3. DESCRIPTION OF THE OBSERVATIONS 

In both datasets described in this paper (GO- 12499, PI 
Lennon; and GO-12939, PI Sabbi) the Wide Field Channel 
(WFC) of ACS is used in parallel with either the UVIS or 
the IR channels of WFC3 to maximize the efficiency of the 
observations. Since there are no filters in common between 
the WFC3/IR channel and the WFC3/UVIS or ACS/WFC, in 
this paper we will simply refer to the instrument (ACS and 
WFC3), without specifying the channel. 

3.1. The GO- 12499 Dataset 

The //5T program GO-12499 (PI. Lennon) represents the 
first epoch of a proper-motion survey of massive stars in the 
30 Doradus region. The survey covers a projected area of 
^ 14' X 12', corresponding to ^ 210 x 180 pc at the distance 
oftheLMC. 

A total of 15 pointings were used to map the whole region. 
All the observations were taken between 3 October and 29 
October 2011, using WFC3 as primary instrument and ACS 
in parallel. The orientation angle of the mosaic was chosen 
to inc lude the very massive runaway VFTS#16 (lEvans et al.l 
I20T0I) . 

Since the goal of the program is to determine the positions 
of the stars with exquisite accuracy, all the observations were 
acquired in a single filter. Filter F775W was chosen in an 
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Figure 3. HST color-composite image of R136, F775W is in blue, Fl lOW in green and F160W in red. 




Figure 4. Map of tlie deptli of tlie mosaic obtained by stacking all the deep 
exposures of the GO- 12499 dataset. The darker the image, the larger the 
number of overlapping images. WFC3 is up and ACS is down. The image is 
rotated by 90° compared to Figure[T]. 

effort to minimize the effect of the strong and variable nebular 
emission that pervades the entire region. 

Each pointing was observed for one orbit of HST tim, and 
each orbit consisted of one short and four deeper exposures 
with each camera. The exposure times were dictated by the 
buffer dumping of WFC3 and by the length of the orbit. To 
give the survey as much spatial uniformity as possible, obser- 
vations were stepped in such a way that no star would fall in 
the inter-chip gap in more than one of the WFC3 exposures. 

The full dataset consists of 15 exposures of 35 sec, 14 of 
500 sec and 45 of 700 sec with UVIS, and 15 of 32 sec, 
14 of 377 sec and 45 of 640 sec with ACS for a total of 
148 exposures. Once stacked together the deep exposures 
result in a mosaic of ~ 28,000 x 22,000 pixels, with a 



40.00 mas pixel ^^ scale. Figure [T] shows the mosaic, with 
colors added based on ground based data. Figure |4]shows the 
depth map of the deep stacked image. 

The CCD detectors on board HST are subject to a flux of 
energetic particles that progressively damage the silicon lat- 
tice of the detectors and create charge traps that redistribute 
the flux from one pixel to the other during the detector read- 
out process. As a result of this cumulative radiation damage, 
the charge transfer efficiency (CTE) of the CCDs on board 
HST is progressively degradin g. We used the algorithm de- 
scribed bv Anderson & Bedin (i2010l) to correct the effects of 
the degrading CTE directly on the images. This algorithm is 
automatically applied by the ACS calibration pipeline. The 
routine we used to correct WFC3 data for CTE can be down- 
loaded from the WFC3 websiteH- 

3.2. The GO-12939 Dataset 

As was done for program GO- 12499, we are observing the 
field with WFC3 as primary instrument and ACS in parallel. 
Visits are organized in four blocks of visits. Figure |5]shows 
the 36 NIR and 24 UVIS WFC3 pointings superimposed on 
the DSS image of 30 Doradus, together with the ACS paral- 
lels. Each pointing corresponds to a single orbit of HST time 
and consists of a four-step dither-pattern. 

The first and third blocks are observed with WFC3/IR and 
ACS. During each orbit we acquire two exposures in the 
WFC3 Fl low and 2 in the F160W filters, while ACS is used 
to collect four observations in the F658N filter (In Figure |5] 
upper panels, WFC3 is in green and ACS is in magenta). 

-" http://www.stsci.edu/hst/wfc3/tools/cte_tools 
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Figure 5. Upper Panels: Positions of the HST pointings in the filters Fl lOW, F160W (WFC3) and F658N (ACS), superimposed on the 33' X 21' DSS image of 
30 Doradus. WFC3 is in green and ACS in magenta. Lower panels: Positions of the HST WFC3 (F275W and F336W filters) and ACS (filter F555W) pointings. 
WFC3 is in blue and ACS is in magenta. Figure[T]is rotated by 45° with the respect to these plots. 



In the second and fourth blocks, each orbit begins with the 
acquisition of a short exposure in the F336W fiher with WFC3 
and a short ACS image in the F555W fiher We then collect 
two images in the F336W and two in the F275W filters. At 
the same time, we acquire four exposures in the F555W filter 
with ACS (in Figure |5] lower panels, WFC3 is in blue and 
ACS is in magenta). 

Half of the IR observations were acquired between 12 and 
19 December 2012 and are described in the next section. The 
remaining half is scheduled for the summer of 2013. The 
NUV observations will be collected in April and in Septem- 
ber 2013. This observing strategy is designed to exploit the 
natural rotation of the telescope during the year 

The filters, the number of exposures, and the exposure times 
used in HTTP are summarized i n Tabl jT] Updates on the sta- 
tus of the survey can be found at |http://30dor.stsci.edu 

Observations with very low background (< 10 e^ for 
WFC3 and < 50 e~ for ACS, such as those acquired in the 
NUV and/or through narrow-band filters) suffer large losses 
for very faint sources, to the point that faint sources can com- 
pletely vanish during the readout transfers. Since Cycle 20, 
it has been possible to mitigate the effects of the degrading 
CTE by adding an extra amount of photons through an in- 
ternal illumination. This process, called "post-flash", has the 
effect of artificially increasing the background level. All the 
NUV observations of HTTP are acquired using "post-flash". 
We decided to not use post-flash for the Ha exposures, since 
we estimated that the background would be sufficiently high. 
Even with the background increased, all sources still suffer 
some CTE loss, and it is necessary to make a photometric 
correction for the detected sources. For this reason we will 
correct all the data acquired in the NUV and visual regimes 
using the Anderson & Bedin algorithm. 

4. DATA ANALYSIS 

The analysis of the F775W (GO- 12499) and the NIR 
data was carried out directly on the bias-subtracted and flat- 
fielded exposures processed by the standard HST calibration 
pipeUnes CALACS and CALWF3. The outputs from the HST 
calibration pipelines are still in the raw-detector pixel frame 
and thus retain several kinds of spatial distortions that have to 
be taken into account. 

4.1. Preparation of the Reference Frame 



Observations in the F775W filters cover a projected area of 
^ 14' X 12'. The first step in reducing the data was to create a 
distortion-free reference frame and relate the astrometry and 
photometry of each exposure to this frame. To achieve this, in 
each F775 W deep exposure we measured positions and fluxes 
for all the sources that had more than 100 counts in each ex- 
posure and no brighter neighbors within a 5 pixels radius us- 
ing the program img2xym_WFC . 9x10 (jAnderson & Kind 
20060 This program uses libraries of empirical PSFs to ac- 
count for the spatial variations caused by the optics of the 
telescope and the variable charge diffusion in the CCD and 
creates a geometric-distortion corrected list of sources. 

We used the Two Micr on AH Sky Survey (2MASS) cata- 
log (ISkrutskie et al.ll2()06l) as an initial astrometric reference 
frame. We then found common stars between each image list 
and 2MASS to define the general six-parameter linear trans- 
formation between the distortion-corrected frame of each ex- 
posure into the 2MASS-defined frame. We used these trans- 
formations to collate the individual star lists in the reference 
frame and determine the average position for ^ 110,000 
bright, isolated, unsaturated stars that in the deep exposures 
have three or more coincident detections. 

For the ACS data we initially used the geom etric distortion 
and PSF library of I Anderson & King! dSSS)- Both the PSF 
and distortion solution were determined using data acquired 
before the last Hubble servicing mission (SM4). A careful 
inspection of the position residuals showed that the PSF has 
developed a 2% asymmetry and the distortion solution has 
changed by ~' 4% of a pixel since SM4. Although this is 
a very small change, it is sufficient to introduce noticeable 
errors in very wide mosaics, therefore we used all the bright 
and isolated stars found in the ACS field of view to derive a 
new PSF and update the geometric distortion. Since all the 
observations in this program were acquired with very similar 
rotation angle, we were not able to independently solve for the 
ACS skew term, however we used the common stars between 
ACS and WFC3 to calibrate it relative to WFC3. 

Finally we improved the internal quality of the reference 
frame by iteration. The rm.s. residuals of the average posi- 
tions were less than 0.01 pixel in each coordinate. In this pro- 
cess we have mapped the ACS pixel scale to the size of the 

^' The program was originally designed for ACSAVFC data, but versions 
for both the WFC3 channels are now available. 
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Instrument 


Filter Name Number of 


Exposures Exposure Time 


Post-Flash 


WFC3/UVIS 


F275W 
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YES 


WFC3/UVIS 


F275W 


I 697 


YES 


WFC3/UVIS 


F336W 


I 14 


YES 


WFC3/UVIS 


F336W 2 700 
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ACSAVFC 
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I 13 


NO 


ACSAVFC 
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NO 


ACSAVFC 


F555W : 


5 640 


NO 


ACSAVFC 


F658N 
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NO 


ACSAVFC 


F658N ; 
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NO 


ACSAVFC 


F775W 


I 32 


NO 


ACSAVFC 


F775W 


I 377 


NO 


ACSAVFC 


F775W ' 


I 640 


NO 


WFC3/UVIS 


F775W 


I 35 


NO 


WFC3/UVIS 


F775W 


I 507 


NO 


WFC3/UVIS 


F775W ' 


I 690 


NO 


WFC3/IR 


FllOW 


I 300 


- 


WFC3/IR 
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I 640 


- 



WFC3/IR 
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I 
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Figure 6. Stars found in a 360 x 360 pixel region in the F775W filter. Sharp 
filaments and unresolved knots of gas are visible. Different colors and sizes 
of the circles correspond to different finding iterations in the photometry. 

WFC3 UVIS channel, so that in all the catalogs a pixel cor- 
responds to 0.040". The final reference list was then used to 
compute the final astrometric transformation of each exposure 
(for both ACS and WFC3) into the reference frame. 

4.2. The Source-Finding Routine 



The stars used in Section HTTI to align the single exposures 
to the reference frame were found in a single pass through 
each image. This approach will find almost all the bright stars 
in a field, but will miss many of the obvious faint sources, 
especially those in the wings of bright objects. 

To identify the fainter sources and push the photometry as 
deep as p ossible, we used an ev olution of the program de- 
scribed in [Anderson et al.l (120081) . The new program will be 
fully described in a forthcoming paper (J. Anderson et al., in 
preparation). 

The previous version of the code was designed to work only 
with ACS and could handle only two filters. This new routine 
can handle multiple instruments and several filters at the same 
time. The routine first finds the bright stars, then subtracts 
them and searches for fainters stars in the subtracted images. 



It performs several such finding iterations (Figure |6), using 
detection thresholds that depend on the number of exposures 
available at each point in the field. This approach is particu- 
larly useful when, as in the case of our survey, the coverage is 
not uniform (in some regions we have only one or two expo- 
sures, while in others we have up to 14 overlapping exposures, 
FigurelHi. An inspection by eye shows that the code efficiently 
avoids spurious detections caused by the diffraction spikes of 
saturated stars and sharp gaseous filaments (Figure |6). 

We were able to identify more than 1.1 x 10^ sources in the 
F775W observations. The photometry was calibrated to the 
VEGAmag photometric system using the Zero Points listed 
on STScI webpageQ These Zero Points are derived for im- 
ages combined using "astrodrizzle", thus they cannot be di- 
rectly applied to our fit-based photometric catalogs. To cal- 
ibrate the fit-based catalogs into the Vegamag system we 
selected several isolated bright stars in the drizzled images 
and measured their magnitudes using 0.5"aperture photome- 
try for ACS and 0.4"for WFC3 data and applied the appropri- 
ate Zero Point. For these stars we then determined the differ- 
ence between the fit-based magnitude and the magnitude in 
the Vegamag system and used the average differences as the 
final Zero Points. 

The F775W photometric catalog was used to measure the 
sources in the IR images. Of the initial 1.1 x 10® sources, 
^ 480, 000 are in the region observed with the first block of 
IR data. 

In addition to finding and measuring stars, the reduction 
routine also produces several diagnostics to help identify the 
stars that are well-measured. One such parameter is Q, which 
reports the linear-correlation coefficient between the PSF and 
the stellar image: Q= 1.0 is a perfect fit. This parameter is 
particularly useful in identifying unresolved objects such as 
background galaxies and blended sources. We also used the 
rm.s. of the individual photometric measurements about the 
mean to help identify well measured stars. Figure |7]shows Q 
as a function of magnitude for the two F775W filters (ACS 
and WFC3 are showed separately) and for the NIR data. 

In both ACS and WFC3 images acquired with the F775W 
filters, stars brighter than mpTTsw ~ 15 are saturated even in 
the shorter exposures. As a result the Q parameter for these 



^^ http://www.stsci.edu/hst/wfc3/phot_zp_lbn 

|http://www.stsci.edu/hst/acs/analysis/zeropoints 
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Figure 7. From the top to the bottom, clockwise, Q parameters as a function 
of magnitude for the filters F775W (ACS and WFC3 respectively), F160W 
and FllOW. Only 10% of all the sources detected in the F775W filers and 
20% of the sources found in the Fl 10 and F160W filters are shown. Objects 
that meet the selection criteria are in black. 

Stars differs significantly from 1. We used aperture photome- 
try to recover the flux of these sources, and we did not apply 
any selection criteria. For ACS we then selected the sources 
in the magnitude range between 15 < mpyvswAcs < 21 with 
Q> 0.95, while for sources fainter than ?71f775Wacs — 21 we 
selected only those objects with Q> 0.92 or photometric error 
< 0.5. Similarly, for the F775WWFC3 filter in the magnitude 
range between 15 < m-F775WwFC3 < 22 we selected only the 
objects with Q> 0.95 and for fainter sources the photometric 
error had to be better then 0.5. 

Because of the coarser resolution of the IR channel, crowd- 
ing is more severe in the filters FllOW and F160W and the 
parameter Q does not provide an estimate of the quality of 
the photometry that is as good as in the case of the optical 
wavelengths. Moreover for most of the field in the NIR we 
have only two measurements per filter, therefore to estimate 
the photometric error we used the signal to noise instead of the 
rm.s. of the individual photometric measurements. As aresult 
for both the Fl lOW and F160W filters we selected only those 
objects with Q> 0.85 and photometric error smaller than 0.1. 
After these selections the combined ACS+WFC3 catalog in 
the F775W filter contains more than 366,000 stars, of which 
- 150, 000 that are detected in both the Fl lOW and F160W 
filters. 

5. COLOR-MAGNITUDE AND COLOR-COLOR DIAGRAMS 

Figure|8]shows the Hess CMDs and CCD for different com- 
binations of the F775W, Fl lOW, and F160W filters. The main 
evolutionary features are highlighted for reference in the top 
left-panel. In this CMD stars brighter than rrapiiow ^ 20 
and bluer than mpyrsw — vtifiiow ~ 0.5 are in the the up- 
per main sequence (UMS). These are intermediate- and high- 
mass stars. The UMS is representative of a stellar population 
younger than ^ 1 Gyr. 

Stars fainter than mpnow ^ 22 and bluer than mpyrsw ^ 
JTiFiiow ~ 1 are in the lower main sequence (LMS). This 
is a mixture of stars with ages spanning the age range from 
potentially very young to several Gyr old. 

Stars brighter than mpnow ^ 21 and redder than 
WF775W — "ipiiow ~ 0.5 are in the sub giant branch (SGB), 



red giant branch (RGB) and red clump (RC) evolutionary 
phases. These sources are characteristic of an evolved old 
(> 1 - 2 Gyr) stellar population. In aU the CMDs, the RC is 
elongated, extending over more than ^1.1 magnitudes in the 
Fl low filter Similarly both the UMS and the RGB are broad 
suggesting that the reddening is variable and differential ex- 
tinction often is significant across the field. 

The lower (mpiiow > 21) right (mpTysw - mpnow > 
1) corner of the CMD is usually populated by young, 
intermediate- and low-mass PMS stars. These sources have 
not yet started the hydrogen burning in their core and are still 
contracting toward the Zero Age Main Sequence. 

In other young (< 5 Myr old) clusters in the MCs, the PMS 
stars are normally well separated from the LMS stars in the 
CMDs derived from HST observations. This is not the case 
for the CMDs shown in Figure |8]not only because differential 
reddening shifts LMS stars toward redder colors, but also be- 
cause with the aging of the stellar population PMS stars move 
towards bluer colors. As the presence of the ^20 — 25 Myr 
old cluster Hodge 301 indicates, star formation in the Taran- 
tula Nebula region has been active for several Myr, therefore 
in the CMDs of the entire field of view, LMS and PMS are 
overlapping. 

5.1. Spatial distribution of the stellar populations 

To better understand the spatial distribution of the identi- 
fied stellar populations, and at the same time, highlight how 
the reddening changes as a function of position, we have di- 
vided our catalog into 18 regions of 3, 000 x 3, 000 pixel each, 
corresponding to '-^ 29pc x 29 pc. Figure |9] shows the loca- 
tion of these regions on the F160W mosaic (Upper Panel) and 
the corresponding 18 mpnow vs mpnow — wpieow CMDs 
(Lower Panel). 

The inspection of these CMDs shows that the UMS is al- 
ways well populated, with the exception of regions A, B, D, 
G, Q, and R. These 6 regions are likely sampling the field of 
the LMC. 

The majority of R136 is in region I, and the NE clump iden- 
tified by Sabbi et al (2012) is in region K. The CMDs for both 
these regions show a large number of red and faint PMS stars, 
as well as very well populated UMSs. Both R136 and the NE 
clump are very young systems (< 5 Myr), and the PMS stars 
are well separated from the LMS stars. PMS stars are likely 
also present in regions L and N. 

Hodge 301 is in region O. The UMS of the corresponding 
CMD is well populated in agreement with the young age of 
the cluster. At these ages PMS stars are too blue to be dis- 
tinguished from the LMS stars using broad-band photometry 
only. 

The RC and the RGB vary from relatively compact and well 
defined (as in regions B, E, M, O, Q, and R) to extremely 
broad. For example in region J, the RC covers more than 
4 magnitudes in the FllOW filter, a typical signature of dif- 
ferential reddening. I n their analysis of the reddening in the 
LMC, bsed on the RC. lHaschke et al.1 ( 1201 Ih find that not only 
30 Dor exhibits the highest reddening of any region in this 
galaxy, but also the highest differential reddening. 

It is interesting to note that young and old stars are affected 
by different amounts of extinction, with the reddening being 
on average less severe for the young stars, than the old ones. 
This effect is particularly evident for example in region H, 
where the UMS is narrow and well defined, but the magni- 
tude of the RC ranges between 18 < mpupw 5. 20. Thi s 
is apparently at variance with conclusions bv Zaritskvl (1 1 9991) . 
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Figure 8. From the top to the bottom, clockwise, CMDs in the filters mpnow vs mpTTSw — "iFllow> "^Fiiow ^s mF775w — "^Fieow. "ipiiow vs 
mpiiow — 'Tipieow. and CCD mpyysw — "^Fiiow ^s /npnow ~ "ipieow- In each diagram both the total number and the maximum number of stars per 
bin of magnitude and color ai'e reported. The most important evolutionary sequences ai'e highlighted in the upper-left panel for reference. 

tinguished by the dust associated to the Tarantula Nebula. We 
will be able to better characterize this issue when we observe 
the region in bluer filters (F275W F336W and F555W). 

The CMDs in Figure |9] suggest that stellar populations of 
different ages are not uniformly distributed. To better high- 
light this effect, in Figure [TO] we have compared the spatial 
distribution of stars selected from different evolutionary se- 
quences. 

In panel A we have selected UMS stars brighter than 
"^Fiiow < 16.6, and bluer than jtihow ^ ^igow < 0.3. 
These stars are younger than ~ 20 Myr and are mostly con- 
centrated around R136 and the NE clump. 

Panel B shows the distribution of stars in the magnitude 
range between 17.9 < rnpiiow < 19.9, and bluer than 
JTiFiiow ~ "ipieow < 0.4. At the distance of the LMC this 
correspond to a mass range between ^ 3 and ^ 7M0. After 
^ 2 Myr the most massive PMS stars in this mass range are 
already merging into the UMS, while a '^ 3 Mq star will still 
be on the MS after '^150 Myr. Therefore panel B shows the 



who found that on average in the LMC old stars are less ob- 
scured by dust than hot young stars^ 

In deriving the reddening map IZaritskvl (119991) used only 
those stars that were detected at the same time in the 4 bands 
U, B, V and I. This selection criterion penalizes the most ex- 
tinguished evolved stars, therefore the average amount of red- 
dening derived for t he RGB stars may h ave been underesti- 
mated. As shown in lHaschke et al.l (12011b in these highly ex- 
tinguished regions the distribution of the reddening values is 
not Gaussian, but shows an extended tail towards the higher 
values. 

It has to be noted that most of the UMS stars found in the 
Tarantula Nebula region are associated with the starburst that 
originated R136 and can be all considered to be at the same 
distance from us. The RGB stars, on the other hand, are 
evenly distributed in the entire body of the LMC and are more 
affected by the 3D structure of the galaxy. If for example 30 
Dor is on the nearer side of the LMC disk, the majority of the 
RGB stars in our field would be beyond 30 Dor and will be ex- 
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Figure 9. Upper Panel Mosaic of the images acquired in the F160W filter. The mosaic has been divided into 18 regions of ~ 3000 X 3000 pixels. The 
corresponding 18 mpnow vs mpnow ~ "ipigow CMDs are shown in the Lower Panel. 

distribution of stars in the age range between ^2 — 150 Myr. 
R136 and the NE clump are still very well defined. The clus- 
ter to the North-West of the image is Hodge 301. A careful 
inspection of panel A shows that Hodge 301 is also visible in 
that map. 

Sources in panel C have colors between 0.7 < mpnow — 
'TiFieow < 1-3 and magnitudes between 20.5 < mpnow < 
22.5. This is the locus of < 1 M© PMS stars. These sources 
are likely younger than ^ 5 Myr. Hodge 301 is not visible in 
this region, while a "chain" of young clustered objects con- 



nects the NE clump to the northern corner of the map. Some 
of these clusters coincide with the embedded massive O stars 
found in the dense knots of dust at the int erface between 
R136 and Hodge 301, already identified bv iBrandner et al.l 
( 120011) . Others are aligned with one side of the shell of ion- 
ized gasJha^un;oimds_a_softand_brigl^^ 
ble dMeaburn II 19841: IWang & Helfandlll991l: iTownslev et alj 
120061) . Although it is possible that the feedback from the two 
massive clusters NGC 2070 and Hodge 301 may have trig- 
gered the star formation in the nearer clumps, the majority 
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of the clusterings are at a projected distance from R136 be- 
tween 30 and 60 pc, too far away to be affected by R136 in 
the less then ^ 2 Myr. Alternatively these systems may result 
from the collapse of residual pockets of gas formed during the 
fragmentation of the giant molecular cloud that formed R136 
(i.e. Bonnell et al. 201 U and references therein). 

LMS stars (panel D) have been selected in the color range 
0.26 < TOpiiow — "ipieow < 0.65 and magnitudes between 
21.5 < mpiiovv < 22.5. These sources have ages between 
a few tens of Myr and several Gyr The spatial distribution 
of these sources is quite uniform in agreement with our as- 
sumption that we are looking at the stars in the field of the 
LMC. The only visible cluster in panel D is Hodge 301, in 
agreement with the fact that old PMS stars cannot be sepa- 
rated from LMS stars using broad-band photometry only. 

The number of LMS stars is significantly lower in the re- 
gions that in panel C correspond to sites of very recent star 
formation. This could be a further indication that 30 Dor is 
on the nearer side of the LMC disk. The spatial distribution 
of RGB stars (panel E, 0.6 < mpnow — "^Fieow < 1-1. 
15.3 < mpiiovv < 18.5) is noisy because of low number 
statistics, but on average appears quite uniform. 

6. DISCUSSION AND CONCLUSIONS 

We have discussed the observing strategy and presented 
preliminary results from the first half of the IR observations 
of the "Hubble Tarantula Treasury Project (HTTP: Unravel- 
ing Tarantula's Web)", an ongoing panchromatic HST survey 
designed to resolve and characterize the stellar populations 
in the Tarantula Nebula down to the sub-solar mass regime 
(< 0.5 M0). HTTP is collecting deep observations of the 
Tarantula Nebula in the NUV (F257W and F336W), optical 
(F555W and F658N), and NIR (FllOW and F160W). These 
data are combined with deep archival observations of the re- 
gion in the F775W filter 

The analysis of the data acquired this far indicates that the 
reddening is highly variable across the region, with the young 
UMS stars being on average less extinguished by dust than 
the old RGB stars. This could be a sign that 30 Dor is on 
the nearer side of the LMC disk. Because the extension of 
30 Dor is likely negligible with respect to the LMC depth, 
the majority of the evolved stars could be beyond 30 Dor, and 
therefore their luminosity would be extinguished by the dust 
associated with the Tarantula Nebula. 

We find that stellar populations of different ages have very 
different spatial distributions. While the intermediate- and 
old-age stars have had time to diffuse in the disk of the LMC, 
the younger stars are still associated with their birth sites. In 
particular the majority of the UMS stars are associated with 
R136, the NE clump and Hodge 301. 

An inspection of the spatial distribution of the PMS stars 
reveals the presence of several small clustered systems. The 
majority of these systems are located between R136 and the 
northern corner of the area covered by our survey, along the 
ridge of an X-ray emitting bubble. We speculate that these 
systems may be minor episodes of star formation that oc- 
curred during the fragmentation of the giant molecular cloud 
that formed R136. Our preliminary analysis confirms that the 
region has actively formed stars for ~ 20 and possibly more 
Myr Once the NUV and optical data are available, we will 
use a Bayesian SED-fitting algorithm to derive the star-by- 
star reddening correction and constrain the age of the various 
sites of recent star formation. A comparison of our data with 
synthetic CMDs will allow us to better assess how star forma- 



tion occurred and is propagating in the Tarantula Nebula, and 
at which star formation rate. 

The uniform coverage and the broad selection of filters of 
HTTP will benefit the community, for example: 

1 . Ha and near IR images will make it possible to iden- 
tify emission-line stars such as Be stars, RGB and SGB 
stars with stellar winds, Herbig Ae/Be stars and low 
mass PMS stars. 

2. Super star clusters such as R136 are often considered 
present-day counterparts of forming globular clusters. 
HTTP offers an observational snapshot of the early evo- 
lutionary phases of these systems and new information 
on how they relate to their surroundings during early 
evolutionary phases. 

3. HTTP will yield the richest and most homogenous 
sample of moderately metal-poor PMS stars of differ- 
ent ages and masses. SEDs and mass-accretion rates, 
which can be derived from this dataset, will provide 
new and much needed constraints on the evolutionary 
models of PMS stars. 

4. HTTP can be used to measure the scale at which UV 
radiation from hot massive stars affects the evolution of 
low-mass accreting PMS stars. 

HTTP is a treasury program and all the data are immedi- 
ately available to the public and can be downloaded from the 
MAST archive. In addition we will deliver a series of high- 
level products to the community, including: 

1 . A unified star catalog for all the filters (F275 W, F336W, 
F555W, F658N, F775W, FllOW, and F160W). For 
each filter we will provide an average flux and an error 
in flux, based on multiple independent observations. 

2. Artificial-star tests done with stars inserted along the 
fiducial sequences. This will allow users to locally 
quantify the level of completeness of the catalog and 
the extent to which crowding may have broadened the 
sequences. 

3. Co-registered stacked images of the field in all filters, 
with calibrated WCS headers. 

4. Maps of differential reddening and the star formation 
history of the field. 

5. A catalog with the properties of all the clusters and as- 
sociations. 

At the moment, only few large ground-based telescopes 
(6m+) can surpass the sharpness of Hubble's data at certain 
wavelengths and on a very limited field of view. In the near fu- 
ture ALMA and JWST will have spatial resolution and sensi- 
tivity similar to HST, but because of the different wavelength 
coverage, they will probe complementary, but distinct stel- 
lar populations and evolutionary stages. In particular JWST 
and ALMA will be able to identify the most dust-embedded 
young stellar objects, highlighting where the most massive 
stars are now forming. In summary, once combined with the 
our ongoing investigations (i.e. stellar dynamics, PI Lennon; 
gas dynamics, PI Gallagher) and studies at other wavelengths, 
HTTP will provide a complete high-resolution picture of the 
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Figure 10. Spatial distribution of the massive UMS stars (panel A), intermediate mass stars in the UMS (panel B), PMS stars (panel C), LMS stars (panel D) 
and RGB stars (panel E). The regions of the mpnow vs mpnow ~ "ipieow CMD used to select the stellar populations are shown in panel F. 



complex interplay between stars, gas, dust and stellar feed- 
back in a bursting regime. HTTP will become the definitive 
catalog of the field and it will provide a unique way to cross- 
identify objects in multiple studies, serving as a touchstone 
for all future works on 30 Dor in particular, and on starbursts 
in general. 
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